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Abstract 

We discuss the "exhaustion" problem in the context of the Periodic Anderson 
Model using Iterated Perturbation theory within the Dynamical Mean Field 
Theory. We find that, despite its limitations, IPT captures the exhaustion 
physics, which manifests itself as a dramatic, strongly energy dependent, sup- 
pression of the effective hybridization of the self consistent Anderson impurity 
problem. As a consequence, low energy scales in the lattice case are strongly 
suppressed compared to the "Kondo scale" in the single impurity picture. The 
IPT results are in qualitative agreement with recent Quantum Monte Carlo 
results for the same problem. 

PACS. 71.28 - Narrow-band systems, heavy fermion metals: intermediate- 
valence solids. 

PACS. 72.15Q - Scattering mechanisms and Kondo effect 
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I. INTRODUCTION 



Metallic compounds containing rare earth elements with partially filled / shells, such 
as CeBe 13 or UPt3, belong to the general category of heavy Fermion materials [|IJ. They 
are characterized by large Pauli susceptibility and specific heat coefficient as compared 
to ordinary metals, which indicate a huge effective electronic mass; and also by anomalous 
transport properties such as non-monotonic temperature dependence of the resistivity These 
anomalies are usually attributed to the formation of a narrow resonant state at the Fermi 
energy due to the admixture of strongly correlated local / orbitals with the metallic band 
of the host. Such materials are normally modelled in terms of an asymmetric Periodic 
Anderson Model(PAM). But many of the high temperature properties of these materials 
are surprisingly similar to those of a Single Impurity Anderson Model(SIAM) |J; so the 
formation of this resonance has also been sometimes interpreted in terms of the SIAM [[|. 
However, an impurity treatment clearly neglects the effects of correlations between the 
impurity sites and also of lattice coherence (i.e. Bloch's theorem) which leads to what 
Nozieres |||| termed as "exhaustion". 

The "Exhaustion" problem |4]||, originally posed by Nozieres in the context of a Kondo 
lattice , occurs when a few mobile electrons, N scr , have to coherently screen many local 
moments, Nf, in a metallic environment, i.e, N scr << Nf. This situation is engendered by 
the fact that only the electrons roughly within Tk (the single impurity Kondo temperature) 
of the Fermi surface can effectively participate in screening the local moments. Thus the 
number of screening electrons can be estimated as N scr = p^e^T^, where pd(e) is the 
conduction band density of states(DOS) and ep is the Fermi level. A measure of exhaustion 
is the dimensionless ratio Q 

N scr Pd(e F )T K 

Nozieres has argued that p is roughly the number of scattering events between a local 
moment and a mobile electron necessary for the mobile electron's spin to precess around 
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a local moment by 2n ||. In the case p > 1, when the number of screening electrons is 
much smaller than the number of local moments to screen, magnetic screening is necessarily 
collective and the single impurity picture becomes invalid. 

In a recent study || of the PAM using the Dynamical Mean Field Theory (DMFT), which 
is exact in the oo— dimensional limit |7|, it was argued that exhaustion was responsible for 
the severe reduction of the Kondo scale of the PAM from that of the impurity problem with 
the same parameters. The issue was explored in detail || using the Quantum Monte Carlo- 
Maximum Entropy Method(QMC-MEM) for solving the self-consistent Anderson Impurity 
problem of the DMFT. The reduction of the Kondo scale and the crossover between the 
two scales was interpreted in terms of Nozieres exhaustion principle ||], and his mapping 
of the model which has exhausted all available conduction band screening states onto an 
effective Hubbard model. Using qualitative and semi-quantitative arguments, Nozieres § 
has suggested that the temperature scale T c associated with the onset of Fermi liquid co- 
herence in this "exhaustion" limit is suppressed compared to the Kondo temperature of the 
single-impurity Anderson Model (SIAM), Tk, by the factor p 

T c ^T K /p = T K p d (e F )/N f . (2) 

While the aforementioned suppression of T c was clearly evident in the QMC simulations, no 
universal relation of T c to Tk such as Eq. ^| was found. 

In this work, we explore the problem of exhaustion in the context of the asymmetric 
PAM once again within the DMFT, but using Iterated Perturbation theory (IPT) || for 
solving the self-consistent impurity problem . We calculate the spectral functions and as 
a check on the accuracy of IPT, compare them with earlier QMC-MEM results ||. We 
find reasonable agreement between the two. The full width at half maximum(FWHM) of 
the "Kondo Resonance" in the spectral functions calculated for the PAM and the SIAM 
provides one measure of the "coherence temperature" T c and the "Kondo Temperature" Tk 
respectively. We find that, consistent with the "Protracted screening" picture of ref |J and 
the "Collective Kondo screening" picture of Nozieres || , a key feature of exhaustion, namely 
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the suppression of T c as compared to Tk, is recovered here as well. Furthermore, we find T c 
and T K to be related as 

T 

Tc = a(U,V)p • (3) 
Here a(U, V) is the fitting parameter and po is defined (similarly to Eq [T]) as po = 
Nf/(p(x)T K ), where p(e) is the bare conduction DOS and x is determined by = 

II. MODEL AND FORMALISM 

The PAM Hamiltonian on a D-dimensional hypercubic lattice is 

+vJ2( d lu + ff.c) + ^ 5> /iT n /iA • ( 4 ) 

In Eq. f|, dj^ (//P) is a creation operator that creates a d(f) — electron with spin a on site i; 
dia{fia) is the corresponding destruction operator. The hopping is restricted to the nearest 
neighbors and scaled as t = t* /2\/T>. U is the on-site Coulomb repulsion for the localized 
/ states, V is the hybridization between d and / states and e/, e d are the site energies for / 
and ti electrons. 

We work in the D — > oo limit where it was shown by Metzner and Vollhardt 0] that 
the irreducible self-energy and the vertex functions become purely local. As a consequence, 
the interacting lattice model can be mapped onto a local correlated impurity coupled to an 
effective bath that is self-consistently determined P,|10||. In this infinite dimensional limit 



the non-interacting DOS has the Gaussian form 0] p(e) = exp 
energy scale such that t* — 1. 
The local / and rf-propagators for the PAM are given by 
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We choose our 



G d , loc (uj) = D L+ - e d - J (5) 
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and G ftloc (u) 
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where D(z) = de^- is the Hilbert Transform of p(e) and a (a;) = u + — ej — T,(u + ) where 

L) + = UJ + zO. 

The self-consistent host is determined by the bare local propagator of the effective single- 
site problem as 

<T 1 = G- } } oc + S (7) 

and self-consistency is achieved through the use of Eqs. |j, ^ and the result for the self-energy 
for the single-site, or the impurity problem £ = £((?) . 

Given a starting self-energy, we use Eqs. |||| and |7| to compute Q, and then a prescription 
for the effective single-impurity problem to calculate the new self-energy. This procedure is 
repeated until self-consistency is achieved. We use Iterated Perturbation Theory(IPT) || 
as the prescription to calculate the self-energy for the effective impurity problem. The 
motivation for using this scheme is that it is semi-analytical and much easier to implement 
than Quantum Monte Carlo(QMC). While it has the disadvantage that it is perturbative,its 
advantages are that we obtain real-frequency data directly at zero temperature and can 
study a wide range of parameters with much less effort compared to the QMC method. We 
now briefly review this scheme. 

The IPT ansatz || for the total self-energy is given by 

A _ U<n f > AE 2 (u) 

where £ 2 (u;) is the second order self-energy defined in terms of a modified bare local prop- 
agator Q^ 1 = Q~ l + €f + /io. The parameter /x is adjusted so as to satisfy the Luttinger's 
theorem @, ImJ^ ^^G f (cu)du = . A and B are chosen so as to reproduce the atomic 
limit and the high frequency behavior of the self-energy at any filling, which yields |§ 

n f (2-n f ) _ 4((l-n f /2)U + e f + no) 

n^-n^Y n (2-n )U 2 ' 1 ' 

Here rif = 2 /"^ pf(u)du is the f-band filling and n is defined by n = 2 /"^ pg(u)duj : 
where pf and are the spectral functions of Gf and Q respectively: Pf = —^lmGf(uj + ), 
and pg = -±Im(j(u + ). 



The conduction band filling is varied by varying e^. We maintain an /-band filling close 
to one, i.e rif ~ 1 by adjusting 6f. The actual value used for rif was 0.99, and the accuracy 
achieved in fixing this value was ~ 1 in 10 5 . For reasons of numerical convenience, we 
calculate the second-order self-energy directly for real frequencies in two steps. We first 
calculate the imaginary part of S 2 using convolution integrals on a lorentzian frequency 
grid; then use Kramer's-Kronig relations to find its real part. Typically we achieve self- 
consistency of the Green functions within 3 to 4 iterations and the solution for ej and /^o is 
found within 10(outer loop) iterations using a non-linear equation solving routine. 

III. RESULTS 

We have performed the calculations described above for various values of U and V and 
we present and discuss some representative results. Since IPT is a perturbative technique, 
we have checked its accuracy by comparing it with the earlier QMC-MEM || results for 
U = 1.5 and V = 0.6. 

Fig. [I] shows the /-spectral function for IPT at zero temperature and QMC-MEM at 
T = 0.025 for three conduction band fillings namely rid = 0.4, 0.6 and 0.8. Both the QMC- 
MEM and IPT results share some common features. For small rid a single narrow Kondo 
peak is centered at the Fermi energy u = 0; however as n<j — > 1, the peak broadens and 
splits in two, with both peaks shifted from the Fermi energy. Apparently, the latter is a 
remnant of the insulating gap found when rif + rid = 2. We see that the IPT results match 
rather well with QMC-MEM for rid = 0.8, but as rid is decreased the deviation increases. 
The difference in the width of the Kondo resonance between IPT and QMC in the low rid 
or the exhaustion limit could be due to two factors: (l)the calculations for the latter were 
carried out at higher temperatures which would lead to temperature broadening, (2) IPT is 
perturbative in U while the impurity Kondo energy scale, Tk is exponential in U. Thus, 
IPT could become less accurate for very small rid since the effective T(0) is decreasing as rid 
decreases (see below). 
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In the QMC-MEM work ||, the SIAM Kondo scales and the PAM coherence scales were 
identified from the T — > limit of the inverse of the appropriate impurity spin susceptibility, 
X7mp(T 0). Since the IPT is not a conserving approximation, it does not provide a unique 
prescription for calculating Ximp- Hence, as an alternate measure of the Kondo and coherence 
scales, we calculate the Full Width at Half Maximum(FWHM) of the Kondo resonance in the 
/-spectral function for the SIAM and PAM and identify these with Tk and T c respectively^. 
We present our results for these in Fig. |](a) for U = 1.5 and V = 0.6. T c is seen to be much 
suppressed as compared to Tk, consistent with Nozieres' arguments HH and the QMC- 
MEM results ||. As mentioned earlier, we find that T c and Tk are related according to 
Eq. [| We have checked this for four sets of parameters and we show the values of a and 
cx(standard deviation of the fit in Eq. [|) i n Fig- @(b) for each of these sets. 

The PAM, even at zero temperature, clearly has many other scales, eg. /t _1 = 
d 2 Yii m /duj 2 \ UJ =o , although they are also related to one other. In the exhaustion limit, 
they are all typically suppressed compared to T K , though to different degrees. We have 
studied all of them, and hope to discuss them elsewhere. 

But the important point is that, given that the PAM is being treated using DMFT, all the 
energy scales in the problem are clearly determined by U and by the effective hybridization 
of the self consistent Anderson-impurity problem Tpam(u) = Im^ 1 = Im [Gjj oc + Sj . 

Tpam is shown in Fig. [|as calculated within QMC-MEM § and IPT, and again the agree- 
ment is reasonable over a fairly wide frequency scale. In as much as Q, the (self-consistent) 
host propagator, has built into it the self-energy, (and hence the "Kondo" screening) pro- 
cesses at all other sites of the lattice except the one under consideration, Tpam(^) nec- 
essarily encapsulates the essential constraints of lattice coherence and collective screening 
that leads to "exhaustion" . In Fig. 01(a), we compare Yp^M^) with the bare hybridization 



1 For the SIAM, it is well known [11] that both Xim P {T — > 0) and the FWHM are proportional to 
the same Tk- 
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Fsiam(u) = 7rV 2 p(u+x)( where p and x are defined below Eq. Note that Tpam is overall 
strongly suppressed compared to T S iam, and in addition has a sharp dip near lu — 0. The 
suppression of T c (and of other scales) compared to Tk is directly related to the suppression 
of Tpam compared to Tsiam , as dramatically brought out by the direct comparison of their 
spectral functions in Fig. £|(b). 

In Fig. |(a), we also compare Tpam with the Hartree hybridization r# artree (u;)(= lm(Gj] oc ) 
where the latter is computed within the Hartree approximation) . r Hartree can be interpreted 
as an effective hybridization that arises when the self-consistent impurity problem is treated 
within the Hartree approximation. Note that it is also strongly suppressed compared to 
T S iam, which can be thought of as being due to exhaustion in the context of "resonant 
level screening" . However, it is rather flat near the fermi level and misses the strong energy 
dependence contained in the dip in T P am(u). As shown in Fig. |](b), the FWHM of the 
Kondo resonance is of course much smaller than that of the Hartree spectral function, 
since the latter has only the resonant-level scales arising from r Hartree, while the former 
includes(at least some of the) non-trivial correlation effects which lead to local moment 
formation and Kondo screening and the associated narrowing of scales. This suggests that 
the strong energy dependent suppression of Tpam is the signature of exhaustion in the context 
of lattice coherent, collective, Kondo screening. 

From the arguments we have presented, one would expect that there should be no such 
exhaustion effects if N s cr >> Nf. We have studied such well, and that this is 

indeed what we find. As we decrease Nf keeping fixed, we find an increase in the FWHM 
of the resonance at the femi level, the dip in the effective hybridization near the Fermi 
level gradually disapears, T Hartree and Tpam, although still somewhat suppressed compared 
to Tsiam, become almost identical and have a weak energy dependence near the Fermi 
level, indicating that the screening processes in this limit are predominantly "resonant level 
screening" . 
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IV. CONCLUSIONS 



We have studied the "exhaustion" problem arising from lattice coherent, collective, Kondo 
screening in the context of the asymmetric PAM in the limit of infinite dimensions within 
IPT. The IPT calculations, despite their limitations, capture( qualitatively and in some 
regimes, even quantitatively) many of the features shown by the QMC calculations ||, 
including the "exhaustion" physics, which manifests itself as a strongly energy dependent 
suppression of the effective hybridization T(u). As a consequence, all the (low) energy scales 
of the PAM are strongly suppressed compared to the corresponding SIAM scales. We have 
presented detailed results for one of them, namely the FWHM of the Kondo resonance, 
and find that its ' suppression ratio is roughly as proposed by Nozieres ( except for a U, V 
dependent scale factor). Together, these results also suggest that IPT would be a useful 
method for incorporating moderate correlation effects into ab initio calculations of Heavy 
Fermion materials [|K| . 
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M. Hettler, D. Logan and Th. Pruschke. This work was supported by NSF grants DMR- 
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S. V and H. R. K). 
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FIG. 1. Comparison of the / density of states obtained from IPT and QMC for three fillings 
of the conduction band, n d = 0.4,0.6,0.8, with n f ~ 1, for U = 1.5, V = 0.6. T = for IPT and 
T = 0.025 for QMC. The agreement is quite good for d-fillings close to one. 
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FIG. 2. (a) The full width at half maximum for the Kondo resonance in the /-spectral function 
, (which is taken as a measure of "T K " for the SI AM and "T c " for the PAM ) for U = 1.5 and 
V = 0.6 as a function of the conduction band filling. The solid line shows Tk/oi(U,V)pq. a is 
the standard deviation for fitting the above relation to T c with a as the adjustable parameter. 
(b)Table.I showing the values of a and a for four sets of U and V. 
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FIG. 3. The effective hybridization within (a) IPT and (b) QMC-MEM for the same parame- 
ters as Fig. |l] showing the decrease in the number of states available for Kondo screening near the 
Fermi level. 
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FIG. 4. (a)The PAM-IPT ,the PAM-Hartree and the SIAM hybridization showing the effects 
of lattice coherence and exhaustion for U = 1.5,V = 0.5, and = 0.4. Near the Fermi level, the 
PAM-IPT r has a dip reflecting a strong energy dependence. Both the PAM-IPT and Hartree 
r are seen to be much suppressed compared to the SIAM T . (b)The corresponding /-spectral 
functions are shown for the same parameters as (a). The Kondo peak broadens for the three cases 
progressively. 
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